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Abstract. We present a fast phase gate scheme that is experimentally achievable 
and has an operation time more than two orders of magnitude faster than current 
experimental schemes for low numbers of pulses. The gate time improves with the 
number of pulses following an inverse power law. Unlike implemented schemes which 
excite precise motional sidebands, thus limiting the gate timescale, our scheme excites 
multiple motional states using discrete ultra-fast pulses. We use beam-splitters to 
divide pulses into smaller components to overcome limitations due to the finite laser 
pulse repetition rate. This provides gate times faster than proposed theoretical schemes 
when we optimise a practical setup. 
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1. Introduction 

Practical quantum information processing (QIP) requires a number of essential 
components outlined by DiVincenzo in pQ. Ion trap schemes fulfil most of these criteria, 
being a particularly promising QIP platform O El [4j [5] . The qubits are represented by 
the internal states of the ions, which permit high fidelity initial state preparation and 
readout j6j Ej. The states of the qubits are manipulated using atom-laser interactions, 
which allow the implementation of a universal set of gates [8j: single-qubit rotations 
are straightforward, and non-separable multiple-qubit gates are mediated by the shared 
motional modes of the ions. The ion trap system has demonstrated its potential with the 
implementation of teleportation [21 EH HI], entanglement of up to 14 ions | [T2l [T3l f!4l fT5j . 
quantum algorithms [EH El CES, EH] and simulation of quantum systems [20l [21] . 

One important criterion is the requirement of a high ratio between the system's 
decoherence time and the timescale for quantum gates. In ion traps decoherence times 
of up to 10 minutes [22] have been achieved, while gates have been performed with speeds 
up to the order of 10/is [22 [24], and fidelity up to 0.993(1) [25]. Most implemented gates 
[SHIEIEHIEHI^ such as the Cirac-Zoller (CZ) [2] and M0lmer-S0rensen (MS) 

[30] gates, excite specific motional sidebands. This restricts the laser intensity that 
can be used and, consequently, limits the gate time to be much slower than the trap 
oscillation period [3T] . For longer, more complex operations, faster and higher fidelity 
gates are required. 

Fast gate schemes have been proposed where all of the phonon modes are excited 
[SU E2J ESI El] , allowing for higher laser power and faster interactions, but complicating 
the disentanglement of motional modes and internal states. More sophisticated 
control is required, using amplitude shaped pulses or carefully timed sequences of 
pulses. Producing the necessary continuous, shaped pulses is a significant experimental 
challenge and such schemes are yet to be implemented. Despite being closer to 
implementation, discrete pulse schemes also face certain limitations. Lasers practical 
for quantum information processing have repetition rates up to around 300 MHz [35] 
and the pulses are fixed in time according to the repetition rate of the ultra-fast laser. 
This is a restrictive, or sometimes prohibitive, limitation for pulsed schemes. 

In this paper we outline a speed-up of quantum gates for ion traps through 
pulse sequences that are readily realised experimentally. We consider a pulse splitting 
technique using unequal path-length interferometers to increase the rate of pulses 
incident on the ions, for a higher effective repetition rate. The pulse splitting process 
also provides certain timing freedoms for the pulse components. These two features 
overcome the hitherto limiting factors for pulsed fast gate schemes. 

The paper is organised as follows: we first consider the fast gate scheme in [33] , 
presented by Garcia- Ripoll, Zoller and Cirac (GZC scheme), which addresses two ions 
and is arbitrarily fast for infinitesimal pulse lengths and separations. The scheme and the 
limitations in its implementation are presented in section [2l In section El we introduce 
the pulse splitting technique and show how it can produce the pulse timings required to 



Fast gates for ion traps by splitting laser pulses 



3 



experimentally implement the GZC scheme. In section HJ we optimise the pulse splitting 
technique and propose a faster gate scheme than existing theoretical schemes for typical 
repetition rates, with a simpler experimental setup. The physical limitations of our 
scheme are outlined in section El 



2. GZC scheme 

We review the proposed GZC fast gate scheme using discrete pulses given in [33]. This 
scheme offers arbitrarily fast gates given perfect laser control, it does not require ground 
state cooling and works outside the Lamb-Dicke regime [36]. The gates are presented 
in the form of a control problem, where two conditions are imposed: 

• A two-qubit phase gate must be completed, which requires the relationship between 
the initial and final qubit states to be described by the gate unitary 

Uj = e iG(T ^. (1) 

We consider the case = 7r/4. 

• The gate should be independent of the motional state, such that ground state 
cooling is not required. For a two qubit system, there are two motional modes: 
the centre of mass and stretch modes. The initial states of both modes (with free 
evolution) should be restored at the end of the gate. 

The solution space for this general control problem is intractable; it grows 
exponentially with the number of pulses. The GZC scheme considers a subset of this 
space, and we now consider the details of the scheme. 



2.1. Gate details 

The Hamiltonian for two ions interacting with resonant light is 

Hj = ^ [a+e ikxi + a+e^ 2 + H.c] , (2) 

where Q(t) is the Rabi frequency, X\ and x 2 are the positions of the two ions, and o\ 
and a 2 are the raising operators for the internal states of the ions. 

We consider TV pairs of counter-propagating n pulses that address both ions and are 
on resonance with the internal transition frequency. Each pulse pair has negligible delay 
between pulses and imparts a 2hk momentum kick on the ions [37], which changes their 
motional state while maintaining their internal states. These simplifying assumptions 
give the evolution operators: 

y[ — JjN ^ e -i2z k rj c (al+a%)(a c +ai) e -iv c 5t k aia c 

_ -QiV -iz k ri r ((Tf-(T%)(a r +al) e -iv r 6t k ala r (3b) 

The subscripts c and r refer to the centre of mass and stretch (or relative motion) modes 
respectively, such that U c is the unitary operator for the centre of mass mode, and a r 
is the stretch mode annihilation operator. The Lamb-Dicke parameter 77 determines 
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r] c = rj/\/2 and r\ r = ry^/4/3, while the trap frequency v = v C) and v r = y/3v. \z\\ is 
the number of simultaneous pairs of pulses incident on the ions at a time Stk after the 
previous group of \zk-i\ pulse pairs, where 5t\ = 0. The direction of the first pulse for 
each pair in the k th group of pairs is determined by sgn (z^). 

Note that the components of the unitaries correspond to rotation and displacement 
operators. For an initial coherent state, then, the motion of the centre of the coherent 
state can be described by a classical trajectory in (x,p)-phase space for each mode. When 
phase space is rotated at the relevant mode frequency to cancel the natural rotation of 
a coherent state, the motional mode restoration conditions require that the phase space 
trajectories are closed. This closed trajectory, determined by the coherent displacement 
operators, maps out an area corresponding to a geometric phase factor. 

Accounting for a global phase factor, our phase gate requires that the relative phase 
shift between the computational excited and ground states (|00), |11)) and the partly 
excited states (|01), 1 10) ) is tt/2. The unitaries in equation f l3aj) show that the excited 
and ground states are only displaced in phase space for the centre of mass mode, while 
the partly excited states are only displaced in stretch mode phase space. We thus require 
that the relative difference between the geometric phase factors for each mode is given 
by 7r/2 to satisfy our control problem. 

Using the evolution operators, the control conditions can be quantified: 

N m-1 

© = 4?7 2 ^2 Z mZ k 

m=2 k=l 

N 

C c = J>,e-^=0 (4b) 
k=i 

N 

a = 5Z^e-^=0, (4c) 
k=i 

where the time of incidence of the k th group of pairs is and 5t m k = t m — tk > 0. 
Equation ( Hal) represents the difference between geometric phase factors as discussed. 

The solution space still grows exponentially with the number of pulses, thus the 
schemes satisfying the condition equations presented in [33] include simple symmetries. 
We refer to the GZC gate scheme as the solution to the condition equations f |4aj) 
presented in [33] which offers arbitrarily fast gate times given arbitrarily fast laser pulses 
on demand. The scheme is characterised by instantaneous groups of pulse pairs z sent 
at times t, interspersed with free evolution: 

z (— 2n, 3n, — 2n, 2n, — 3n, 2n) (5a) 
t = (-7i, -r 2 , -r 3 , r 3 , r 2 , n). (5b) 

As above, at time tk € t, Zk G z determines the number of simultaneous counter- 
propagating pulse pairs \zk\ that are incident on the ions, as well as the direction sgn(z^) 
of the first pulse in each pair. Every pulse pair in a group has the same first pulse 
direction. To obtain arbitrarily fast gate times, the number of pulses in the scheme, 



sm(is5t mk ) 



sin(y / 3^t mfe ) 



7T 



(4a) 



Fast gates for ion traps by splitting laser pulses 
N = 14n, is scaled towards infinity. 



5 




Rotating Xc Rotating Xc 



(a) n=4 (b) n=32 

Figure 1: Phase space trajectories for the GZC gate, where Xc and Pc are the 
dimensionless position and momentum for the centre of mass mode. Here phase space is 
rotating with frequency v. Increasing numbers of pulses, n, reduce the time evolution, 
seen through smaller angles in the second plot. The larger kicks are also apparent, 
maintaining the required geometric phase. The initial state and trajectory direction are 
shown. 

Figure [U shows the rotating phase space trajectory plots for the centre of mass 
mode with increasing n. Phase space trajectories identify the behaviour of a given 
solution scheme, as well as providing a check on the scheme validity since the closure 
of trajectories and geometric phase are apparent. Note that as the number of pulses, 
or kicks, is scaled towards infinity, the angles (evolution times between kicks) shrink as 
the kick sizes increase, preserving the area and maintaining the closed trajectory. The 
scheme's star shape aids in this effective scaling with n. 

2.2. GZC limitations 

Gate times are restricted by dissipation in the system, due to coupling to the 
environment, as well as control errors, caused by factors such as noise in the pulse 
amplitude and duration, a non-harmonic trap and evolution between the 'instantaneous' 
pairs of pulses. The most significant factor, however, is the experimental limitation on 
the repetition rates of the necessary ultra-fast lasers. 

As n is scaled up, we require increasing numbers of pulses to be incident on the ions 
in a timescale much smaller than the trap period. A typical maximum repetition rate 
for lasers useful for QIP is given by 300MHz [35]. This restricts the number of pulses 
where the small timescale approximation holds and the error for the scheme is less than 
10 -4 , a significant threshold for error correction [3 [25]. At this repetition rate, the GZC 
scheme as presented in [33] does not achieve the requisite fidelity even for n — 1, while 
an alternative fast gate scheme [32j is limited to gate times on the order of the trap 
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period 284ns [27] . A faster repetition rate of 1GHz permits the GZC scheme for n = 1 
at twice the trap period, however the n = 2 case still has too large an error from the 
assumed instantaneous pulses. 

A further restriction is imposed on the timing of ultra-fast laser pulses. The laser 
emits pulses at regular time intervals, while the GZC scheme requires the pulse group 
times tk to be free. In the next section we will present a pulse splitting technique, 
making the scheme compatible with the timing and repetition rate restrictions. 

3. Splitting scheme 

Recently, Senko et al. used pulse trains, where large pulses are split into smaller 
components, to apply a fast spin-dependent momentum transfer to an atom [38]. This 
motional control is performed in the strong excitation regime, where the Rabi frequency 
Q is greater than the trap frequency v. The pulse splitting turns each laser-emitted 
pulse into multiple ion-incident pulses using beam-splitters, overcoming the repetition 
rate restriction. The technique also gives some pulse timing freedoms; the path lengths 
for the different pulse components after splitting can be varied. This technique can be 
used for significant gains in implement able gate times. 

There is an immediate limit on the pulse splitting technique: we are restricted by 
the maximum pulse area (largest Bloch sphere rotation) that we can generate with our 
laser. If each laser pulse can contain energy sufficient for up to 1024 n pulses, using the 
splitting technique this corresponds to a maximum of 512 pairs of tt pulses. This is still 
a significant advantage: the effective repetition rate (incident on the ions) is 512 times 
the repetition rate of two counter-propagating lasers without pulse splitting. 

3.1. Application to GZC scheme 

The GZC gates become experimentally viable (although complex) using the pulse 
splitting technique, since both the repetition rate and timing restrictions are relaxed. We 
present a possible optical setup for the GZC scheme in Figure El With pulse splitting, 
the quasi-instantaneous pulse groups can be produced by time delays much smaller than 
the repetition period using beam-splitter loops, as shown in the figure. The number of 
pulses in the GZC scheme scale with n. The figure represents n — 1, while n = 2 J can 
be implemented by introducing extra beam-splitter loops identical to F, connected to 
the front of the optical setup (directly before the (4/7): (3/7) beam-splitter) as Y is to 
the T\ — r 3 delay loop. We will assume that the delay time for such a loop is zero so that 
the pulse groups arrive simultaneously. In practice there will be a small delay producing 
negligible error (details in section [5]). For loops labelled with delays, such as T\ — r 3 , the 
marked delay corresponds to the total delay as compared to the straight optical path. 

This complex setup would be experimentally challenging to implement, since it was 
not designed with the pulse splitting technique in mind. It is possible to use alternate but 
similarly complicated beam-splitter and mirror arrangements to perform the scheme. In 
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Figure 2: Optical setup for the GZC gate, with n = 1. The pulse delay times, permitted 
by tailoring path lengths, are marked on the delay loops such that the GZC scheme is 
performed. Loops without explicit delays, such as Y, are assumed to have zero time 
delay. The Roman numerals on particular pulse paths corresponds to the ordering of the 
pulses. In the path through A, each pulse is split in two at X. Both of the component 
pulses are subsequently split in two, with delays such that pulse components I and II 
arrive as a pair, then III and IV (with I arriving immediately before II, and III before 
IV). The path through B involves a complex setup before point Z where the delays 
are as short as possible such that a pulse is split into thirds: three component pulses 
of equal area reach point Z. These pulses are split into ordered pairs as with the path 
through A. The beam-splitters are 50/50 unless otherwise marked with fractions. The 
different paths permit the pulse group numbers, timings and directions as specified by 
the GZC scheme. 
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the following section we remove the GZC symmetries and consider schemes designed for 
experimental simplicity with the pulse splitting technique, optimising the pulse timing 
and direction freedoms. 

4. Optimal schemes 

4-1- Cost function 

To find optimal solutions to the condition equations, we must first introduce an 
appropriate cost function. Our cost function J rewards fast gate times Tq ) balanced 
with high fidelity, as quantified in section 15.11 We write our cost function in terms of 
the error E — 1 — Fp, where F P is the process fidelity [39] : 



where A and B are constants introduced to provide a balance between the gate time 
and fidelity in the cost function. Simulations using various values of A and B indicated 
that appropriate values are A — 10, B = 100. The solution space for this cost function 
is not convex, as shown in Figure^ thus finding the global minimum is a hard problem. 
The figure shows the log of the above cost function for clearer contours, and the true 
cost function has the same features. A useful function involving the condition equations 
with a convex solution space could not be determined. The figure shows a 2D-slice of 
the solution space for the GZC scheme with n — 1, where r 3 is set to its solution value 
associated with the global minimum. The third dimension further obscures the global 
minimum, as do extra dimensions for more free variables in later explored schemes, thus 
searching for solutions (optimal gate schemes) becomes increasingly computationally 
difficult. 

Using our cost function, J, we optimised the gate times for different experimental 
setups using the NLopt nonlinear-optimization package [40], implementing a controlled 
random search with local mutation, as outlined in [4T1 142] . 

4.2. Simple scheme: direct pulses 

We consider first the simplest setup for the pulse splitting technique. Ideally, we want 
to use a single laser to avoid mismatches in pulse timing and frequency, introducing 
new sources of error. Since we require pairs of counter-propagating pulses, the simplest 
setup is in the form of Figure BJ The pulse pairs are formed at the last beam-splitter, 
X. One of the paths after this beam-splitter is slightly shorter than the other, such 
that the pulses are incident separately on the ions (with minimal separation in time for 
each pair). The first pulse in each pair always arrives from the same direction using 
this scheme, while in the following section we consider alternating the direction of the 
first pulse in each pair. We require resonant tt pulses to satisfy the simplified controlled 
phase gate conditions in equation f l4al ). 




(6) 
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Figure 3: Solution space shape for log (J) for our cost function J, given the GZC scheme 
with n — 1. The r 3 value is set to 1.8T P (trap periods), its solution value. The global 
minimum has been circled. 

Figure [4] illustrates the pulse splitting process for each laser pulse split into eight 
components, or eight splittings, which corresponds to four counter-propagating pulse 
pairs. More splittings are easily produced: the number of splittings is doubled when an 
extra beam-splitter loop, of the form shown with a t\ delay in the figure, is introduced. 
Each added loop links to the others in the same way the t\ delay loop is joined to the 
i 2 delay loop. These delay times are optimised to find the fastest gate schemes. 

The optimal gate schemes for this setup provide gate times on the order of the trap 
period; the fastest gate found is 1.37Tp (389ns), where T P is the trap period. This gate 
is for one pulse split into eight pairs, and the next fastest gates are around twice the 
trap period. We use a Lamb-Dicke parameter of r\ = 0.2 [27], and a repetition rate of 
300MHz. 

It is significant to note that each solution set of optimised delay times for different 
numbers of pulses and splittings contained an exact multiple of 0.5 and of This is 
equivalent to a kick with direction ( — 1) in the rotating centre of mass and stretch mode 
phase space trajectories respectively. It follows that the positive direction of each pulse 
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Figure 4: Optical setup for splitting each pulse into 8 pulses, in 4 pairs, with variable 
delays t\ and t 2 . Each pulse is split into a pair of pulses at X, after which one component 
pulse is incident on the ions from above in the figure, and its counterpart is incident 
from below. The path lengths are tailored such that the second pulse in each pair is 
incident on the ions directly after the first. 



pair kick results in closed phase space trajectories only when the evolution time permits 
an effective kick reversal. This explains the restriction of the scheme timescale to the 
order of the trap period. We thus consider a scheme with little extra complexity that 
permits the direction reversal of the first pulse in a pair. 

4.3. Simple scheme: alternating pulses 

When a direction switch is introduced for certain pairs of pulses, as in Figure El there 
is a small speedup in the gate time to below the trap period. Two splitting loops with 
slightly different lengths (i 2 and t 2 + e in the figure) allow the ordering of the counter- 
propagating pulses to be reversed. Figure [6] shows the gate times for different numbers 
of pulses. Different numbers of splittings of each emitted pulse are compared, where 
each number of splittings corresponds to an experimental setup, with more splittings 
introducing extra beam-splitter loops. 

Note that each setup approaches a gate time lower limit around 0.8T P . Larger 
numbers of splittings, corresponding to more power applied to the ions, approach this 
limit faster. As for the direct pulse scheme, we are restricted by the trap period 
timescale, since each solution involves delay times 0.5 + = 0.79T P . The lower 
limit increases slowly with the number of pulses, as the laser repetition rate is taken 
into account. In principle, the fastest gate corresponds to sending one pulse with very 
large area, which is split into a large number of pulses performing a phase gate. In 
practice, the number of splittings is restricted by the laser capabilities. For a given 
maximum pulse energy, the number of possible splittings is fixed and the fastest gate is 
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Figure 5: Optical setup for splitting each pulse into four pairs of pulses, with variable 
delays t\ and £ 2 - The numbers correspond to the arrival order for the pulses, as for the 
GZC scheme. Note that the first pulse of the first pulse pair arrives from the top, while 
the first pulse of the second pulse pair, number 3, arrives from the bottom. Thus the 
pulse pairs have two possible direction orderings. 

achieved by some optimal number of pulses, determined by the minimum for each curve 
in Figure EH Our optimisation considered up to 128 pairs of pulses, which requires 11 
beam-splitters. 

The limiting timescale indicates that introducing a reverse direction pulse pair on 
a single delay loop doesn't provide the required freedom for schemes with scalable gate 
times as the applied power is increased. In phase space, this represents the requirement 
for this alternating scheme to have trap period delay times to close the phase space 
trajectories. The GZC scheme is not limited by the trap timescale, which means that 
appropriate placement of the reverse pulse pairs removes the timescale limitation. We 
identify optimal reverse pulse pair timings through a free search of solution space, and 
reverse engineer the practical experimental setups from these solutions. 

4-4- Reverse Engineering 

The three condition equations suggest that three variables are required for solutions, 
however added free variables can introduce further flexibility and reduce the cost 
function for more optimal schemes. It was found for the simple setup schemes in the 
previous two sections that even when extra delay loops, and corresponding variables, 
were added, only three of the delay times took nonzero values for optimal solutions. 
This suggests that the symmetry of the scheme plays a significant role. Intuitive design 
of simple optical setups with gate times that scale with power proved challenging. 
We instead use a free search of solution space and identify symmetries leading to 
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Figure 6: Gate operation time for different numbers of pulses, using the alternating pulse 
direction scheme. Each curve corresponds to a different number of splittings of each laser 
pulse, corresponding to the indicated number of split pairs. The curves approach a line 
with a slight gradient, which represents the limiting gate time for the scheme. The slope 
of the line is governed by the laser repetition rate, 300MHz. 



implementable setups. The limitations of global optimisation searches restrict our search 
to low numbers of free variables. 

Free variable searches were carried out for optimised times t = (0, Xi, x 2 , x 3 , £4, x 5 , ...) 
associated with kicks z = (l,— 1,1,— 1,1,— 1,...). A search over five free variables re- 
vealed that the optimal solution takes the following symmetric form with (a, 6, c) = 
(1,1,1): 

t = (-7i, -r 2 , -r 3 , r 3 , r 2 , n) (7a) 
z = (an, — 6n, cn, — cn, 6n, —an). (7b) 

As before, \zk \ is the number of pairs of pulses applied to the ions with first pulse direction 
sgn(z/e) at time This (1, 1, 1) scheme, though straightforward to implement, is still 
limited by the trap period, approaching a gate time around 0.73Tp. 

Nine free variables give more variety in optimal solutions. For 8 split pairs, a 
symmetric scheme in the form of equation f lTaj ) with (a, 6, c) = (1, 2, 2) was found to be 
optimal. Higher numbers of split pairs gave non-symmetric optimal solutions, such as 
a scheme for 32 split pairs (320 total incident pulse pairs) with a gate time of 0.086Tp. 

The non-symmetric schemes are extremely complex to implement using the beam- 
splitter loops. It also becomes hard to find the equivalent schemes for higher numbers 
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of splittings, since the global optimisation search is not deterministic, and we have a 
complicated high-dimensional search space. In contrast, the general symmetric form in 
equation ( 17aj) provides a lower-dimensional search space, and certain values of (a, 6, c) 
for the symmetric form permit straightforward experimental setups. Note that this 
symmetric form includes the GZC solution, which takes the values (a, 6, c) = (2, 3, 2). 

These symmetric schemes are often well suited to solving the condition equations 
such that the gate time is not limited by the trap period, as for the GZC 
case. The schemes, although sometimes suboptimal, are among the better solutions 
for the free variables searched. Searching low integer values for (a, 6, c) to 
find implement able, scalable schemes yielded the (1,2,2) scheme as the optimal 
experimentally implementable scheme. Figure [7] shows the optical setup of the (1, 2, 2) 
scheme. It is similar in structure to the GZC setup, with some complexity removed. 

As derived in [34] , the optimal relationship between applied pulse pairs TV and gate 
time Tq follows the power law Tq oc N~ 2 ^ . The (1, 2, 2) scheme and the GZC scheme 
both follow this scaling for large numbers of pulses, with a slightly better prefactor 
(15% smaller) for the (1, 2, 2) scheme as shown in FigureEl Lower numbers of splittings 
approach this scaling at different rates for either scheme, seen in the figure. 

The scaling suggests infinitesimal gate times as the number of incident pulse pairs 
increases towards infinity. The setup for each scheme involves a single pulse split into the 
required components for the entire gate operation, which means the number of incident 
pulse pairs are directly linked to the laser emitted pulse energy. The prior limitations in 
fidelity due to finite laser repetition rates are thus transferred to a gate time limitation 
due to finite laser powers. 



As well as the power cap, we discuss in section 15.21 that errors in each pulse can 
have a significant impact on the gate fidelity. We have discussed multiple schemes 
with different setup complexities and scaling relationships. The fastest scheme for less 
than 20 incident pulse pairs is the simple alternating pulse scheme (section 14.31) , which 
provides a gate time of 1.18Tp for 16 incident pairs, given a very simple setup. For 
20 or more incident pulse pairs, some extra complexity gives gate times scaling most 
efficiently with the laser power for the (1, 2, 2) scheme. This scheme achieves a gate time 
of 0.12Tp (33ns) for 320 incident pulse pairs. This is more than two orders of magnitude 
faster than experimental gates with ion trap systems [23, 124] . 

5. Limitations 

Gate fidelity loss is caused by imperfect laser control, as well as dissipation from coupling 
to the environment. The effects of environmental decoherence are considered in [34J, 
and were found to cause exponential decay of the fidelity. However, the dissipation 
limit on gate speeds was found to not restrict our scheme. We find that control errors 
introduce a stronger bound on gate speed. Imperfect laser control encompasses the 
effects of finite pulse durations, non-instantaneous evolution between pulses, mistimed 
pulses, and errors in the amplitude and duration of each pulse. 
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Figure 7: Optical setup for the (a, 6, c) = (1, 2, 2) scheme for n = 2, i.e. 20 incident pulse 
pairs. The roman numerals mark the order of incidence for pulses travelling through the 
marked path. Beam-splitter loops without a marked time delay represent a negligible 
delay. The loop marked Y doubles the number of splittings; without it we would have 
the n — 1 case, while another such loop would give us n = 4. 



5.1. Pulse timing 



Using the condition equations f j4"al!4bll4cl) , we determine the significance of shifting pulses 
in time, analysing the effect of mistimed pulses and pulses assumed to be instantaneous 
(using zero as a beam-splitter loop delay time). The error from failing to satisfy the 
motional conditions is derived in |34j : 

d-exph^Cel 2 ] (8) 
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Figure 8: Gate times for the (1, 2, 2) symmetric scheme and the GZC scheme for different 
numbers of pulses incident on the ions. Both schemes are performed by a single laser 
emitted pulse split into the required incident pulse pair components. The linear fits 
represent the power law scaling for each scheme, gate time Tq oc N~ 2 ^ 3 for N pulse 
pairs. The scaling coefficients are 6.30 and 5.37 for the GZC and (1,2,2) schemes 
respectively. 



C 2 ^ exp[-^|?7 r CV| 2 ] 



E m = (Q- C\ 



4dC 2 )/8. 



(9) 
(10) 



Here we have assumed a pure initial internal state, separable from an initial low energy 
thermal motional state. Ideally, the internal and motional states become entangled 
during the gate, before the components are again separable upon the gate's completion. 
The error is found using the evolution of the internal state density matrix coefficients, 
detailed in [34J. 

The upper bound phase error for the process fidelity when the acquired phase 
= 7r/4 — x, for some error \x\ < 7r/4, is given by 



E P <~ 4 



cos(2x). 



(ii) 



The error E = E m + E p given here is defined as 1 — F P) where F P is the process 
fidelity. The total error for failing to satisfy the conditions in Eqs. f |4al ) is thus given by 



E < - \ cos(2x)) + \ (6 - Ci - C\ - AC.C,) 



(12) 
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The pulse timings were systematically shifted, and it was found that for an error 
upper bound of 10 -4 , the system was stable for systematic time shifts of beyond 14ps 
(using a repetition rate of 300MHz, trap frequency of 2tt x 3.52MHz). Pulse duration can 
be much smaller than this value, and this corresponds to a systematic mirror placement 
error of 0.4cm, which is much larger than the accuracy achievable in experiments. The 
pulse time shift errors are thus not a limiting factor for this scheme. 



5.2. Pulse area 

Pulse duration and amplitude errors, on the other hand, can be significant. Both 
duration and amplitude errors affect the area of the pulse sent, which affects the area 
of each pulse component in the splitting scheme. The incident pulses on the ions will 
not be exact tt pulses, meaning that the condition equations are no longer applicable. 
We consider the worst case process fidelity F w [39j : 

F w = mm(\(4>\UjU e \i;)\ 2 ), (13) 

where tj) is the quantum state of the ions, thus we are minimising over motional and 
internal states. U e is the realised unitary operator for the process with errors, and Uj is 
the ideal gate scheme unitary operator, 

U I = e i ^ a *. (14) 

Recall that this ideal unitary is composed of momentum kicks for each ideal pair 
of pulses [z\ = ±1): 

U Ip = e -2i**i(*i*f+*2*i) ? (15) 

interspersed by the free evolution of the system, which rotates the motional modes in 
phase space: 

JJ _ e -iv8t k (aia c +aia r ) (^6) 

For a general atom-laser interaction, if 6 is the rotation about the Bloch sphere 
caused by the pulse of duration T, i.e. 6 — J Q ^-dt : the unitary is given by 

U kick = (cos(0) - ism(6)a*e- lkxi(T i)(cos(6) - isin(0)ofe- <fcB2<T 5) (17) 

To examine tt pulses with some small error, we consider 9 = | + e, e<C |. The 
unitary U p for a pair of counter-propagating pulses can be expanded in terms of e: 

U p = U Ip + Y,c j fi, (18) 

3 

where each fj is found from equation f fTTI ). The unitary for a pulse scheme is given by 
combinations of free evolution, [7 r , and pairs of pulses, U p . For a four pulse pair scheme: 

U € = U p U r U p U r U p U r U p . (19) 

We considered the worst case fidelity to second order in e, as the first order terms 
cancel. For the four pulse pair scheme, a lower bound on the worst case fidelity was 
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found to be F w = 1 — 331 e 2 . This means that for an error of 1% in the tt/2 pulse energy, 
for example, e = 5 x 10 -3 and the worst case fidelity is 99.2%. 

The fidelity can scale very poorly as the number of pulses increases, as expected for 
systematic errors in the pulse area. Substituting equation ( IT5I) into equation f [T9l) gives 
four terms (fC) to first order in e: 

U e = U z + e(h{U r U Ip f + U Ip U r h{U r U Ip f + {U Ip U r ) 2 hU r U Ip 

+ (U Ip U r ff 1 ) + 0(e 2 ). 

This corresponds to permutations of fi and Ui p . Thus for an TV pulse scheme there are 
±C = N such terms. Similarly there are TV second order terms of U e corresponding to 
permutations of / 2 , and second order cross terms in f\. Higher orders m of e will 
contribute cross terms in fi. For large AT, this scales as N m . Note that for e on the 
order of 0.01, and TV on the order of 100, perturbation theory breaks down, indicating 
that the gate errors become large. 

Pulse area errors are thus a limiting factor for high fidelity fast gates. We have 
considered an absolute worst case bound, and experimental error is expected to be 
much less than this bound. The scaling of error with pulse numbers is nonetheless 
concerning. These errors are not unique to pulse splitting schemes; the pulse area 
errors must be considered for any gate scheme requiring a given pulse area and multiple 
pulses. Techniques such as stimulated Raman adiabatic passage (STIRAP) [2H [43] 
exist for limiting pulse control errors, which would greatly improve the fidelity of gate 
schemes, although possibly at the cost of gate speed. Even with low pulse numbers, we 
can achieve gate times faster than the trap period. 

5.3. Pulse direction 

The angular stability of the lasers also plays a part in the fidelity of the scheme. 
Systematic errors were introduced independently for the two pulse directions to represent 
misalignment with the axis of the ion trap. As one would expect, the sensitivity to pulse 
direction becomes more significant with the number of pulses. The (1, 2, 2) scheme for 
20 pulses has error less than 10 -4 for angular precision of 0.36°, or 6.3 milliradians. The 
scheme with 80 pulses requires stability to 0.09°, or 1.6 milliradians. Again we see that 
errors are compounded from increasing the number of pulses. 

6. Conclusions 

We have considered gate schemes optimised for both gate time and experimental 
simplicity using the pulse splitting technique. The most practical optical setup schemes 
are limited by the trap period timescale, and achieve gate times of up to 0.8T P (0.23/xs). 
Some added complexity provides gate times that scale with the number of pulses applied 
to the ions. Our (1, 2, 2) scheme is limited by control errors such as laser intensity and 
duration fluctuations, however gate speeds faster than the trap period can be achieved 
with high fidelity. For 80 incident pulse pairs, the scheme achieves gate times of 0.29Tp 
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(82ns), more than two orders of magnitude faster than experiments, and more pulses 
give faster gate times. For comparison, an alternate fast gate scheme [32j, using directly 
incident sequences of n pulses, has an operation time longer than 100ns for the same 
300MHz repetition rate. 

The techniques explored for the gate implementation can be extended to consider 
the excitations of higher numbers of ions by introducing new control parameters to 
account for the added motional modes. Alternatively, the gate could be implemented 
in a shuttling architecture [44]. When thus extended, the gate scheme presented would 
be implementable for more complex quantum information processing operations. 
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